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摘要 :【 目的】 建立 花椒 窒 吉 丁 Agrilus zanthoxylumi j&, % fik f& $h 3 28 24 3E E , 3 AER e T8 X SL ,为 
今后 研究 其 触角 的 化 学 感受 机 制 及 生物 防 控 提 供 理论 支撑 。【 方 法 】 采 用 高 通 量 测 序 平台 Tllumina 
NovaSeq 6000 对 花椒 罕 吉 丁 雌 雄 成 虫 触角 进行 转录 组 测序 ,用 Trinity 软件 对 获得 的 高 质量 reads 进 
行 序列 拼接 与 组 装 ;使 用 BLAST 软件 将 触角 转录 组 数据 比 对 NR, NT, Swiss-Prot, GO, KEGG, 
BLASTX, eggNOG, Pfam, TmHMM, SignalP, KO, Map, BLASTP 和 RNAMMER 公共 数据 库 ; 基 于 初 
步 算 选 到 的 花椒 罕 吉 丁 候选 气味 结合 蛋白 (odorant binding protein, OBP) 和 化 学 感受 蛋白 
(chemosensory protein, CSP) ZU J- 463) A FG k 45 8] 78 4& €) 85 HE BRUT 91,4070. MEGA 软件 进行 
系统 进化 分 析 。 运 用 RPKM (reads per kilobase per million mapped reads) 值 对 嗅觉 相关 基因 表达 量 
进行 分 析 。 【结果 】 花 椒 窒 吉 丁 唆 雄 成 虫 触 角 转 录 组 测序 共 获得 36 209 条 基因 和 90 982 条 转录 
本 ,其 N50 分 别 为 2 103 和 2 523 bp, 组 装 完整 性 较 高 。 注 释 到 NR 数据 库 的 基因 最 多 (41.62% ) , 
其 中 与 赤 拟 谷 盗 Tribolium castaneum 相似 基因 所 占 比 例 最 高 (19% ) 。 在 GO 数据 库 中 比 对 到 
11 614 个 基因 , 按 功 能 分 为 细胞 组 分 、 分 子 功能 与 生物 学 进程 三 大 类 57 个 分 支 ,其 中 分 子 功能 大 类 
中 的 结合 (70.57% ) 与 催化 活性 (45.51% ) 相关 基因 占 比 最 多 ;KEGG 代谢 途径 分 析 表 明 ,7 427 条 
基因 参与 了 5 类 代谢 通路 ,其 中 涉及 信号 转 导 的 基因 最 多 ,为 815 条 ;筛选 到 7 个 候选 OBP 基因 和 
5 个 具有 全 长 开放 阅读 框 的 CSP 基因 ,其 编码 蛋白 均 具 有 化 学 感受 蛋白 家 族 的 典型 特征 。 系 统 进 
化 分 析 表 明 ,花椒 窄 吉 丁 OBPs 和 CSPs 分 别 与 苹果 小 吉 丁 A. mali 的 OBPs 和 CSPs 氨基酸 序 列 一 
致 性 最 高 。RPKM 值 表明 ,嗅觉 相关 基因 AzanOBP1 和 AzanOBP2 在 只 成 下 触角 中 不 表达 ,在 雄 成 
下 触角 中 微量 表达 ;4zanrOBP3 在 雄 成 虫 触 角 中 高 丰 度 表达 。 【结论 】 首 次 获得 了 花椒 窄 吉 丁 成 下 
触角 转录 组 数据 ,筛选 到 了 花椒 窄 吉 丁 OBP, CSP, Or, IR 和 SNMP 等 的 嗅觉 相关 基因 。 推 测 触角 
中 高 丰 度 表达 的 OBPs 对 雄 成 虫 识别 同类 异性 释放 的 信息 素 或 寄主 植物 释放 的 挥发 物 起 关键 作 
用 。 研 究 结果 可 为 花椒 窗 吉 丁 化 学 感受 基因 功能 分 析 及 嗅觉 感受 机 制 研究 黄 定 分 子 基 础 。 

关键 词 : HE SOT; 触角 转录 组 ; 基因 注释 ; 气味 结合 蛋白 ; 化 学 感受 蛋白 
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Analysis of the antennal transcriptome and olfaction-related genes of 


Agrilus zanthoxylumi ( Coleoptera: Buprestidae ) 

GONG Xue-Fang, XIE Shou-An* , YANG Ping, GUO Li, CHEN Di, CHE Xian-Rong, WANG Yan- 
Lai, LÜ Shu-Jie (College of Forestry, Northwest A&F University, Yangling, Shaanxi 712100, China) 
Abstract: [ Aim] To establish the antennal transcriptome database of Agrilus zanthoxylumi adult and to 
explore target genes related to olfaction, so as to provide the theoretical basis for the future study of the 
chemosensory mechanism of its antenna and its biological control. [ Methods] High-throughput 


sequencing platform Illumina NovaSeq 6000 was used for transcriptome sequencing of the antennae of 
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female and male adults of A. zanthoxylumi. "Trinity software was used for sequence splicing and assembly 
of the high-quality reads obtained. BLAST software was used to compare the antennal transcriptome data 
with public databases NR, NT, Swiss-Prot, GO, KEGG, BLASTX, eggNOG, Pfam, TmHMM, SignalP, 
KO, Map, BLASTP and RNAMMER. Phylogenetic analysis was carried out using MEGA software based 
on the nucleotide sequences of the preliminarily screened candidate odorant binding proteins ( OBPs) and 
chemosensory proteins ( CSPs) of A. zanthoxylumi, and homologous proteins from other coleopteran 
insects. The RPKM ( reads per kilobase per million mapped reads) value was used to evaluate the 
expression levels of olfaction-related genes. [Results] By sequencing of the antennal transcriptome of 
female and male adults of A. zanthoxylumi, a total of 36 209 unigenes and 90 982 transcripts were 
obtained with the N50 of 2 103 and 2 523 bp, respectively, indicating the high completeness of 
assembly. Most unigenes (41.6296) were annotated to the NR database, with the highest proportion 
(1996 ) of genes similar to those of Tribolium castaneum. In the GO database, 11 614 unigenes were 
annotated and divided into 57 branches of three functional categories, cellular component, molecular 
function and biological process, among which binding (70. 5796) and catalytic activity (45. 5196) 
related unigenes in the category of molecular function were the most. The KEGG metabolic pathway 
analysis indicated that 7 427 unigenes are involved in five categories of metabolic pathways, with the most 
unigenes (815) involved in signal transduction. Seven candidate OBP genes and 5 CSP genes with full- 
length open reading frame were screened, and their encoded proteins have the typical characteristics of 
chemosensory protein family. The phylogenetic analysis showed that the OBPs and CSPs of A. 
zanthoxylumi have the highest sequence identity with the OBPs and CSPs of A. mali, respectively. The 
RPKM values showed that AzanOBP1 and AzanOBP2 were not expressed in the antennae of female adults, 
but slightly expressed in the antennae of male adults, and AzanOBP3 was expressed with high abundance 
in the antennae of male adults. [ Conclusion] For the first time, the antennal transcriptome data of A. 
zanthoxylumi were obtained, and the olfaction-related OBP, CSP, OR, IR and SNMP genes of A. 
zanthoxylumi were screened. It is speculated that OBPs expressed with high abundance in the antennae 
play a key role in identifying pheromones released by female adults of the same species or the volatiles 
released by host plants. The results can lay a molecular basis for further studying the functions of 
chemosensory genes and olfactory sensitivity mechanism of A. zanthoxylumi. 

Key words: Agrilus zanthoxylumi; antennal transcriptome; gene annotation; odorant binding protein; 


chemosensory protein 





i6 WUTE d; T Agrilus zanthoxylumi Jg 85 38 H 
( Coleoptera) zi J | ( Buprestidae ) "E di ] JB Agrilus, 
是 我 国花 椒 林 中 最 主要 的 毁灭 性 性 干 害 虫 ( 吴 海 ， 
2006; 权 丽 君 ,2016 ) 1987 年 在 陕西 省 韩 城市 椒 
园 首 次 发 现 花 椒 罕 吉 丁 是 致使 花椒 大 量 死 亡 的 主要 
原因 ( 薛 耀 宗 ，1989 ) 。 目 前 , 主要 分 布 于 陕西 、 甘 
肃 ` 山 东山 西 等 地 (党 国 刚 等 , 2017 ) 。 幼 虫 主要 性 
食 枝 干 形成 层 , 引起 树 势 衰弱 ,产量 品质 下 降 ( 李 
J£, 2013) , 主要 危害 特征 是 受害 部 位 虫 疤 松 软 湿 
润 , 有 无 色 或 淡 褐 色 液 体 渗 出 ,俗称 “ 流 油 ”( 焦 浩 ， 
2011)。 成 虫 取 食 叶片 补充 营养 ,造成 叶片 的 缺 刻 、 
空洞 , 树 皮 流 胶 、 软 化 腐烂. 干枯、 龟 裂 最 后 脱落 至 
整 株 死亡 ( 杜 平 , 2012; ERSE, 2012), IE Ti 
























































丁 的 研究 主要 集中 于 虫 情调 查 种群 动 态 、 发 生 规律 
( 李 孟 楼 等 , 1990) 及 其 化 学 防治 ( 康 克 功 等 ,2007 ) 
等 方面 。 近 年 来 ,在 花椒 窄 吉 丁 与 寄主 植物 互 作 的 
分 子 机 制 方面 ,对 寄主 挥发 物 .后 肠 与 疮 便 提 取 物 对 
其 引诱 和 驱 避 作用 进行 了 初步 研究 ( 袁 丽 芳 等 ， 
2016; xJ2 SE, 2017) ; ATHERE T J ME HER E 
fib fg 25 pr de A AED B ETI T jn ye (DU y 55, 
2019) ;对 气味 结合 蛋白 基因 AzanOBP3 的 克隆 、 原 
核 表 达 及 组 织 表 达 谱 分 析 也 进行 了 初探 ( 巩 雪 芳 
等 , 2020) ,但 是 对 其 作用 的 分 子 靶 标 尚 不 清楚 。 
触角 是 昆虫 感受 外 界 信 号 的 主要 器 官 ,在 长 期 
进化 过 程 中 形成 了 能 够 感知 各 种 化 学 信息 物 的 高 度 
灵敏 且 特 异 的 嗅觉 感受 系统 ( 王 桂 荣 等 ,2002; 修 
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伟 明 等 , 2005; Jo SE, 2012; Brito et al., 2016) , 
EME dz SOR BIOS URR BB ARE HE TER 
衍 过 程 中 均 发 挥 重要 作用 (Zhang et al., 2015) , IR 
觉 识别 过 程 需 要 多 种 蛋白 的 参与 ,主要 有 人 气味 结合 
蛋白 (odorant binding proteins, OBPs) 、 化 学 感受 和 蛋 
H (chemosensory proteins, CSPs) 气味 受 体 (odorant 
receptors, ORs) 、 离 子 型 受 体 (ionotropic receptors, 
IRs) .感受 神经 元 膜 和 蛋白 (sensory neuron membrane 
proteins ，SNMPs ) 气味 降解 酶 (odorant-degrading 
enzymes, ODEs) (JE XA, 2018) FE RR AEk 
巴 液 内 发 现 两 类 作为 环境 信息 分 子 结合 载体 的 水 溶 
性 和 蛋白 , 即 气味 结合 蛋白 和 化 学 感受 蛋白 (Angeli et 
al., 1999; Yang et al., 2011), 

花椒 窗 吉 丁 是 如 何 感知 寄主 气味 物质 并 准确 定 
位 到 唯一 的 寄主 植物 ,在 寻找 寄主 过 程 中 嗅觉 是 如 
何 发 挥 作 用 的 ? 深入 研究 花椒 窄 吉 丁 嗅觉 系统 有 助 
于 阐明 其 定位 寄主 及 寻找 配偶 时 识别 信息 物质 及 其 
化 学 通讯 与 感受 机 制 ,而 获得 大 量 嗅觉 相关 基因 是 
进行 分 子 机 制 研 究 的 前 提 。 为 此 ,本 研究 借助 高 通 
量 测序 对 花椒 窄 吉 丁 成 虫 触角 转录 组 进行 生物 信息 
学 分 析 ,获得 了 参与 嗅觉 感受 的 一 些 重要 化 学 感受 
基因 OBP, CSP, OR, IR fil SNMP 家 族 基 因 ,为 今后 
花椒 窄 吉 丁 嗅觉 感受 蛋白 的 作用 机 制 提 供 理论 文 
撑 , 对 开展 害虫 行为 调控 剂 的 筛选 及 绿色 防 控 技术 
研究 具有 重要 的 意义 。 



















































































1 材料 与 方法 


1.1 供 试 昆 虫 

2018 年 5 -7 月 于 陕西 省 蓝田 县 普 化 镇 花椒 园 
采集 花椒 鹤 吉 丁 成 虫 后 带 回 实 验 室 于 25% 下 饲养 。 
选取 1 000 头 活 动能 力 较 强 的 成 虫 ,用 灭 菌 三 蒸 水 
冲洗 虫 体 表 面 ,70% 的 酒精 表面 消毒 30 s, KA= 
水 冲洗 虫 体 2 次 后 解剖 镜 下 于 冰 上 迅速 切 下 触角 肉 
雄 分 开 收集 , 液 氮 冷 冻 保 存 于 -807C 备用。 
1.2 雌雄 成 虫 触角 总 RNA 的 提取 与 质量 检测 

采用 Trizol 法 分 别提 取 花 椒 罕 吉 丁 肉 、 雄 成 忠 
触角 总 RNA。 利 用 1. 096 琼脂 糖 凝 胶 电 泳 和 凯 奥 
K5500 分 光 光 度 计 检 测 RNA ak £ ( OD444/0D», FE 
值 ) ,安捷伦 2100 RNA Nano 6000 Assay Kit 检测 完 
整 性 和 浓度 ,质量 合格 的 RNA 样品 保存 于 - 807€ && 
用 或 者 立即 进行 反 转 录 。 
1.3 文库 构建 与 测序 

1.2 节 样 品 RNA 质量 及 纯度 经 检测 合格 后 ,用 





























带 有 Oligo( dT) 的 磁 珠 富 集 mRNA ,加 入 片段 缓冲 液 
$ mRNA 中 使 其 片段 化 成 为 短片 段 ,以 片段 后 的 
mRNA 为 模板 ,用 六 碱 基 随机 引物 合成 cDNA 第 1 
链 , 并 加 入 缓冲 液 、dNTPs、RNaseH 和 DNA Polymerase 
I 继续 合成 cDNA 第 2 链 , 采 用 QIA Quick PCR iX 
剂 盒 纯 化 并 加 EB 缓冲 液 洗 脱 。 对 洗 脱 纯化 后 得 到 
的 双 链 cDNA 再 进行 末端 修复 .加 碱 基 A、 加 测序 接 
头 处 理 , 回 收 目的 大 小 片段 进行 PCR 扩 增 。 在 获得 
clean reads 后 需 采 用 Trinity 软件 对 clean reads 拼接 
和 组 装 ,获取 分 析 所 需 的 参考 序列 。 转 录 本 具有 可 变 
剪接 的 特点 ,组 装 出 全 长 转录 本 ,对 测序 得 到 的 原始 数 
据 (raw data) 过 滤 , 获 得 高 质量 的 数据 (clean data) ,在 
所 有 的 转录 本 中 , 取 最 长 的 contigs 作为 unigenes。 
1.4 转录 组 数据 分 析 

高 通 量 测序 中 ,Q20 值 是 指 测 序 时 碱 基 识 别 过 
程 中 ,对 所 识别 的 碱 基 给 出 的 错误 概率 ,错误 识别 概 
率 是 1% ,正确 率 99% ;质量 值 Q30 ,错误 识别 概率 
0. 196 ,正确 率 99.9% ;N50 即 履 盖 50% 核 背 酸 的 最 
大 unigene K EREA 9i] 38 i HE J& HE, unigenes 
序列 与 NR, NT, Swiss-Prot, GO, KEGG, BLASTX, 
eggNOG, Pfam, TmHMM, SignalP, KO, Map, 
BLASTP 和 RNAMMER 公共 数据 库 比 对 ,对 获得 的 所 
有 unigenes 进行 功能 注释 及 结果 统计 。 
1.5 嗅觉 相关 基因 的 鉴定 与 序列 分 析 

BLAST 同 源 搜索 从 触角 转录 组 中 鉴定 OBPs, 
CSPs, Ors, IRs 和 SNMPs 嗅觉 相关 基因 。 从 NCBI 
数据 库 中 下 载 苹果 小 吉 丁 4 mali MAIE J A. 
planipennis 等 其 他 鞘翅 目 昆 虫 的 嗅觉 相关 基因 核酸 
序列 ,采用 ClustalX 进行 多 序列 比 对 ,利用 MEGA 软 
件 邻 接 法 (neighbour-joining method, NJ method ) 构 
建 系统 发 育 树 ,对 结果 进行 1 000 次 bootstrap 验证 。 
1.6 嗅觉 相关 基因 的 表达 分 析 

衡量 基因 表达 量 时 ,单纯 以 map 到 的 read 数 计 
算 基因 的 表达 量 , 统 计 上 是 不 合理 的 。 本 研究 利用 
每 百 万 条 reads 中 来 自 于 某 基 因 每 千 碱 基 长 度 的 
reads ZX ( reads per kilobase per million mapped reads, 
RPKM ) ( Mortazavi et al., 2008) ,对 嗅 党 相关 基因 表 
达 量 进行 分 析 。 
2 结果 
2.1 测序 结果 及 从 头 组 装 

花椒 罕 吉 丁 雌 雄 成 虫 触角 转录 组 数据 
(SUB6796283 ) 经 过 滤 后 分 别 获 得 43 124 393 和 
45 854 648 条 clean read ,组 装 得 到 36 209 条 基因 和 



























































1162 昆虫 学 报 Acta. Entomologica Sinica 63 35 





90 982 条 转录 本 ,转录 本 和 基因 的 N50 分 别 为 2523 ”上 ,Q20 H 10096 , Q30 均 不 少 于 9296 ,组 装 完整 
和 2 103 bp (图 1)。 样 品 clean data 均 达 到 9696 以 性 较 高 ( 表 1)。 


3000 
2500 
2000 
1E -£ 1500 
E 
1000 
500 
0 
SS Od d Qd edu d d du 
S S CLEA LEL LLC ENEON AS S SIS 
S a er C SS Na g A rA SY Fa 7 
WOW 9 SS GE Ve e v X Nr S Xe uS P ME AS es 
序列 长 度 (bp) 
Sequence length 


图 1 花椒 窗 吉 本 雌雄 成 虫 触角 转录 组 unigenes 长 度 分 布 


Fig. 1 Length distribution of unigenes in the antennal transcriptome of female and male adults of Agrilus zanthoxylumi 




















X1 dEMUETET BH DL BL fe SUBE AR 注释 结果 做 比较 韦 恩 图 (图 2), 共 有 36 209 条 
Table 1 Summary of assembly for the antennal unigenes 得 到 注释 ,其 中 包括 29 927 条 (54. 2196 ) 
E of female and malg adults of 肉 成 虫 触角 的 和 26 706 4& (50.78% ) 雄 成 虫 触角 
grilus zanthoxylumi : 4 f 
rm MUERE HJ, NR 数据 库 注 释 到 的 unigene 最 多 (15 070 条 ， 
Number Total number of bases MESES 41. 62% ) » 同 时 注 TE 到 NT, NR, BLASTX 和 
转录 本 Transcripts 90 982 66 565 475 39.18 BLASTP 四 大 数据 库 的 有 2 682 条 。 同 JE 白 的 物 
基因 Unigenes 36 209 183 683 635 37.63 . Nn - 
一 一 一 一 一 一 一 一 一 一 种 来 源 分 析 结果 表明 ,与 未 拟 谷 盗 的 相似 基因 最 多 
2.2 基因 功能 注释 (1996) (E13) . 


使 用 BLAST 软件 将 获得 的 花椒 窗 吉 丁 瞧 雄 成 
虫 触角 转录 组 unigenes 序列 与 NR, Swiss-Prot, GO, 
BLASTX 和 BLASTP 等 公共 数据 库 比 对 ( 表 2), 并 对 
表 2 花椒 窄 吉 丁 雌雄 成 虫 触角 转录 组 基因 功能 
注释 结果 统计 
Table 2 Statistics of gene function annotation in the 


antennal transcriptome of female and male adults of 
Agrilus zanthoxylumi 














数据 库 Unigene 数量 di EE C96 ) 
Database Number of unigenes Proportion 
NR 15 070 41.62 ^ F. 
NT 3 692 10.19 saatt 
Sviss-Pro 12 508 34.54 2 TEIR T MEHR MR P AEE 
GO 11 614 32.08 注释 结果 韦 恩 图 
BLASTX 11 676 32.25 Fig. 2 Venn diagram of gene annotation results of 
eggNOG 8 766 24.21 the antennal transcriptome of female and 
Pfam 9 428 26.04 male adults of Agrilus zanthoxylumi 
TmHMM 2 500 6.90 "T 
SignalP 1 301 3.59 2.3 GO 功能 分 类 
M ESEO iSc GO 功能 分 类 统计 (图 4) 表 明 , 比 对 到 11 614 
Map 4 025 11.12 条 unigenes ,分 为 细胞 组 分 (cellular component) 、 分 
BLASTP 9 458 26.12 子 功 能 (molecular function ) 与 生物 学 进程 (biological 


el B 9.o process) 三 大 类 共 57 个 分 支 。 其 中 ,生物 学 进程 大 
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赤 拟 谷 盗 Tribolium castaneum 

旺 犀 金 包 Oryctes borbonicus 

中 欧 山 松 大 小 圳 Dendroctonus ponderosae 
WER Lasius niger 

WEAK EL Zootermopsis nevadensis 

其 他 /未 知 Others/unknown 














多 3 NR 数据 库 中 注释 的 花椒 罕 吉 丁 雌 雄 成 虫 触角 unigenes 的 物种 分 布 


Fig. 3 Species distribution of unigenes in the antennal transcriptome of female and male adults of 


Agrilus zanthoxylumi annotated in the NR database 


类 中 参与 细胞 进程 (cellular process ) 基因 最 多 
(72.55% ) , 其 次 是 代谢 进程 (metabolic process ) 
(58. 1596) , ^E Eg 5€ Wi] 5 ( biological regulation ) 
(42. 9196 ) ,其 余 均 在 40% 以 下 ;细胞 组 分 大 类 中 细 
胞 部 分 (cell part) 基因 最 多 (78. 54% ) ,其 次 是 细胞 
器 (organelle ) ( 49. 19% )、 细 胞 器 部 分 (organelle 
part) (40. 6996 ) 、 膜 (membrane ) (30. 4896 ) 与 膜 部 
分 (membrane part) (27. 00%). 大 分 子 复 合 物 
( macromolecular complex ) (26. 4796 ) , 其 余 均 低 于 
10% ;分 子 功能 大 类 中 结合 (binding ) (70. 5796) 5 
催化 活性 (catalytic activity) (45. 5196 ) 的 基因 最 多 ， 
其 余 均 在 10% 以 下 。 
2.4 KEGG 代谢 通路 分 析 

通过 整合 数据 库 , 分析 unigenes 可 能 参与 或 涉 
及 的 代谢 途径 的 注释 结果 (图 5) 发 现 有 7 427 条 
unigenes 参与 了 5 类 代谢 通路 : 细胞 进程 (cellular 
process) ,环境 信息 处 理 (environmental information 





processing ) , 遗传 信息 处 理 ( genetic information 
processing) , 新陈代谢 (metabolism) 和 有 机 系统 
(organismal system) 。 其 中 unigenes 注释 最 多 的 代 
谢 通路 为 信号 转 导 (signal transduction) (815 条 ) , 翻 
译 (translation ) (586 2& ) , 内分泌 系统 (endocrine 
system) (466 条 ) ,运输 和 分 解 代 谢 (transport and 
catabolism) (440 条 ) , Hr Æ , 2] TARIFE fit ( folding, 
sorting and degradation) (417 条 ) , 碳水 化 合 物 代谢 
(carbohydrate metabolism) (397 条 ) 。 
2.5 嗅觉 相关 基因 的 筛选 及 进化 分 析 

初步 鉴定 到 7 个 OBP 基因 (Azan0BP1 -7) ,其 
中 只 有 5 个 OP 基因 具有 全 长 开放 阅读 框 ， 
AzanOBPl, AzanOBP2, AzanOBP3, AzanOBPA 和 
AzanOBP5 具有 6 个 保守 的 半 胱 氮 酸 位 点 ,符合 气味 








结合 蛋白 的 特征 ( Cys-X15-39-Cys-X3-Cys-X21-44- 
Cys-X7-12-Cys-X8-Cys,X 代表 任意 一 种 氨基 酸 ) ,为 
典型 的 OBP。5 个 具有 全 长 开放 阅读 框 的 CSP 基因 
(AzanCSPl -5) ,3 个 OR 基因 (hzanOR1 -3),5 ^ IR 
基因 (hzanIR1 -5) 和 5 个 SNMP 基因 (AzanSNMP1 - 
5)。 基 于 AzanOBP 和 AzanCSP 基因 以 及 其 他 类 群 
昆虫 OBP 和 CSP 基因 核 昔 酸 序列 构建 系统 发 育 树 ， 
发 现 花椒 罕 吉 丁 的 7 个 OBP 基因 分 支 聚 类 较为 分 
散 ( 图 6) ,但 都 至 少 与 一 个 鞘翅 目 昆 虫 的 OBP 同 源 
基因 聚 类 在 一 支 , 且 序 列 一 致 性 很 高 。4zan0BP1 和 
AzanOBP3 的 序列 一 致 性 较 高 , 且 分 别 与 白蜡 罕 吉 
Agrilus planipennis AplaGOBP83a ( GenBank 登录 号 : 
XM. 025974145. 1) TI SE SE /h Si T Agrilus mali 
AmalOBP2. (GenBank 登录 号 : MF615491. 1) 聚 到 一 
X AzanOBP2 和 AzanOBPA 均 与 苹果 小 吉 丁 的 
AmalOBP3 (GenBank 登录 号 : MF615492. 1) 亲缘 关 
系 近 ,三 者 聚 到 同一 支 。 

花椒 罕 吉 丁 5 个 CSP 基因 聚 类 较为 集中 (图 
7) ;AzanCSPl , AzanCSP2 和 AzanCSP3 聚集 到 同一 
分 支 ,5 个 CSP 基因 均 与 苹果 小 吉 丁 的 AmalCSP 亲 
缘 关 系 最 近 。 对 于 鞘翅 目 其 他 不 同 种 的 昆虫 的 CSP 
基因 ,也 出 现 一 定 的 聚 类 关系 分 别 在 进化 树 中 形成 
多 个 独立 的 小 分 支 ;同一 种 昆虫 的 不 同 CSP R 
类 关系 较 分 散 也 形成 了 多 个 独立 的 小 分 支 。 
2.6 嗅觉 相关 基因 的 表达 量 

利用 RPKM 值 评估 不 同 嗅觉 相关 基因 在 肉 雄 
成 虫 触角 中 的 表达 量 , 结 果 表 明 ( 图 8) , AzanOBPI 
和 AzanOBP2 在 雌 成 虫 触 角 中 不 表达 ,在 雄 成 虫 触 
角 中 微量 表达 ;Azan0BP3 在 雄 成 虫 触角 中 高 丰 度 表 
达 ( RPKM > 1 000) ,在 肉 成 虫 触 角 中 中 量 表达 。 
AzanOBPA4 在 雄 成 虫 触角 中 中 量 表达 (RPKM > 














1164 


昆虫 学 报 Acta Entomologica Sinica 


63 卷 





分 子 功 能 Molecular function 


细胞 组 分 Celluar component 


生物 学 进程 Biological process 


Protein tag 

Nutrient reservoir activity 
Molecular carrier activity 
Hijacked molecular function 
Translation regulator activity 
Antioxidant activity 
Molecular transducer activity 
Signal transducer activity 
Molecular function regulator 
Structural molecule activity 
Transcription regulator activity 
Transporter activity 

Catalytic activity 

Binding 


Cell 

Nucleoid 

Virion 

Virion part 

Other organism part 
Synapse 

Supramolecular complex 
Synapse part 
Membrane-enclosed lumen 
Cell junction 
Extracellular region 
Extracellular region part 
Macromolecular complex 
Membrane part 
Membrane 

Organelle part 

Organelle 

Cell part 


Biological phase 

Cell aggregation 

Detoxification 

Reproduction 

Cell killing 

Pigmentation 

Presynaptic process involved in chemical synaptic transmission 
Rhythmic process 

Cell proliferation 

Growth 

Signaling 

Biological adhesion 

Immune system process 
Behavior 

Locomotion 

Multi-organism process 
Reproductive process 
Multicellular organismal process 
Response to stimulus 
Localization 

Developmental process 

Cellular component organization or biogenesis 
Biological regulation 

Metabolic process 

Cellular process 


0 10 20 30 40 50 60 70 80 90 


基因 百分比 Percentage of genes 


4 TERET T MERE Uff unigenes 的 GO 分 类 


Fig. 4 GO classification of unigenes of female and male antennae of Agrilus zanthoxylumi 


400) , AzanOBP 基因 在 雄 成 虫 触角 中 的 表达 量 排 序 
为 : ManOBP3 (RPKM = 1 160. 57) > AzanOBP4 
( RPKM = 244. 93 ) > AzanOBP6 (RPKM = 145. 09) > 
AzanOBP7( RPKM = 44. 78) > AzanOBP5 ( RPKM = 
42.36) > AzanOBP1 ( RPKM = 2. 30) 和 AzanOBP2 
( RPKM 21.79) EMERE ft ff rS BUHEIY DJ : 
AzanOBP3 ( RPKM = 205. 25) > AzanOBP6 ( RPKM = 


42.94) > AzanOBPA ( RPKM = 39. 69) > AzanOBP7 
( RPKM = 30. 66) > AzanOBP5 ( RPKM = 14. 04) > 
AzanOBP1 ( RPKM -0) fil AzanOBP2 ( RPKM 20) ( [El 
8: A)。 除 AzanCSP1 ( RPKM > 200 ) 在 雄 成 虫 触角 
中 中 量 表达 外 ,其 他 4zanCSP 基因 在 雌雄 成 虫 触 角 
中 表达 量 均 较 低 (RPKM < 200); AzanCSP2 Te WE JV, 
虫 触角 中 不 表达 (RPKM =0) ,在 雄 成 虫 触角 中 微量 
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图 5 HEIT T EHE d fff unigenes 的 KEGG 代谢 通路 


Fig. 5 KEGG metabolic pathways of the antennal unigenes of female and male adults of Agrilus zanthoxylumi 


























A: 细胞 进程 Cellular process; B; 遗传 信息 处 理 Genetic information processing; C; 环境 信息 








谢 Metabolism; E: 有 机 体系 统 Organismal system. 


表达 (RPKM =3.064) ; AzanCSP 基因 在 雄 成 虫 触 角 
中 的 表达 量 排序 为 : AzanCSP1 ( RPKM = 302. 16 > 
AzanCSP5 ( RPKM = 64. 58) > AzanCSP4 ( RPKM = 
32.07) > AzanCSP3 ( RPKM = 8. 64) > AzanCSP2 
( RPKM =3. 06) , AzanCSP 基因 在 肉 成 虫 触角 中 的 
表达 量 排 序 为 : AzanCSP1 ( RPKM = 207. 85) > 
AzanCSP5 ( RPKM = 79. 65) > AzanCSP4 ( RPKM = 
27.21) > AzanCSP3 ( RPKM = 4. 24) > AzanCSP2 
(RPKM 2-0) (8: B)。 气 味 受 体 基 因 AzanOR 基 
因 在 雌 成 虫 触角 中 的 表达 量 排序 为 : AzanORI 
(RPKM = 3. 89) > AzanOR2 ( RPKM = 0.28) > 
AzanOR3 ( RPKM =0. 18) , AzanOR 基因 在 雄 成 虫 触 
角 中 的 表达 量 排序 为 : AzanOR1 ( RPKM -18. 13) > 
AzanOR3 ( RPKM =2.81) > AzanOR2 ( RPKM =1. 33) 
(图 8: C) 。 离 子 型 受 体 基因 AzanIR Ak A E WER E 
触角 中 的 表达 量 排序 为 : AzanIRA ( RPKM =2. 86) > 
AzanIR2( RPKM =1. 89) > AzanIR3 ( RPKM =1.73) > 
AzanlR1( RPKM =0.95) > AzanIR5 ( RPKM =0.14), 在 





\ 处 理 Environmental information processing; D: 新 陈 代 


雄 成 虫 触角 中 的 表达 量 排序 为 : AzanIR2 ( RPKM = 
6. 08) > AzanIRA ( RPKM = 4. 21 ) > AzanIR1 ( RPKM = 
3.70) > Azanl R3 ( RPKM =2. 08) > AzanIR5 ( RPKM = 
0.48) (图 8:; D) 。 神 经 元 感受 膜 蛋白 基因 AzanSNMPS 
在 肉 成 虫 触角 中 不 表达 (RPKM =0) ,在 雄 成 虫 触角 
中 微量 表达 ( RPKM =0. 24) ,AzanSNMP 基因 在 肉 成 
虫 触角 中 的 表达 量 排序 为 : AzanSNMP2 ( RPKM = 
19. 90) > AzanSNMP1 ( RPKM - 0. 42) > Azan$NMP3 
( RPKM = 0. 05) > AzanSNMP4 ( RPKM = 0. 02) > 
AzanSNMPS5 ( RPKM =0) ,在 雄 成 虫 触角 中 的 表达 量 
排序 为 : AzanSNMP2( RPKM = 105. 17) > AzanSNMP4 
(RPKM - 0. 46) > AzanSNMP1 ( RPKM =0.32) > 
AzanSNMP 5( RPKM =0.24) » Azan$NMP3 ( RPKM = 
0.06) (E18: E), 


3 讨论 与 结论 


嗅觉 相关 基因 的 鉴定 可 为 基于 反 向 化 学 生态 学 
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图 6 应 用 邻接 法 基于 核 苷 酸 序 列 构 建 的 花椒 窄 吉本 雌雄 成 虫 触角 转录 组 中 7 个 AzanOBP 基因 














与 其 他 鞘翅 目 昆 贝 OBP 基因 的 系统 发 育 树 (1 000 次 
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Fig. 6 Phylogenetic tree of seven AzanOBP genes in the antennal transcriptome of female and male adults 


of Agrilus zanthoxylumi and OBP genes from other coleopteran insect species based on 


nucleotide sequences by using neighbor-joining method (1 000 replicates) 
基因 来 源 物 种 Origin species of genes; AzanOBP : AE $i ] Agrilus zanthoxylumi ; AmalOBP ; 苹果 小 吉本 Agrilus mali; AplaOBP; HIR T 














Agrilus planipennis; DhelOBP ; 158 2$ P Dastarcus helophoroides; HparOBP ; W SRA fa, Holotrichia parallela; PaenOBP ; Tj 5k E 7E Ut 





























Pyrrhalta 











aenescens; PmacOBP ; if i EUF Pyrrhalta maculicollis; AcorOBP ; 4] 2k WN 4 f& Anomala corpulenta; AchiOBP; EX ^F. Anoplophora chinensis; 
AglaOBP ; 光 肩 星 天 牛 Anoplophora glabripennis; MaltOBP ; 松 黑 天 牛 Monochamus alternatus; TcasOBP ; 赤 拟 谷 盗 Tribolium castaneum; RdomOBP ; 
RR Rhyzopertha dominica; DarmOBP ; 华山 松 大 小 寺 Dendroctonus armandi. 











的 研究 提供 潜在 分 子 靶 标 (Zhu et al., 2017) ,利用 
昆虫 行为 调节 剂 和 分 子 干扰 技术 调控 昆虫 个 体 间 、 
昆虫 与 寄主 间 的 化 学 通讯 来 防 控 害 虫 是 国内 外 研究 
的 热点 ( 杨 海 博 等 , 2018 ) 。RNA-Sed 技术 已 成 为 转 
录 组 研究 的 重要 手段 ,成功 运 用 于 非 模式 昆虫 功能 
基因 挖掘 ( 郑海霞 等 , 2018 ) ,为 深入 研究 昆虫 串 觉 
相关 基因 功能 黄 定 了 坚实 的 基础 。 本 人 研究 对 花椒 罕 






































吉 丁 雌雄 成 虫 触角 转录 组 进行 测序 , 共 获 得 90 982 
条 转录 本 。G0 功能 分 类 表明 ,分 子 功能 分 类 中 参 
与 结合 的 基因 占 比 最 多 (70.57% ) ,猜测 花椒 罕 吉 
丁 成 虫 触角 中 存在 大 量 结合 化 学 信息 物质 的 蛋白 。 
仍 有 59% 的 unigenes 未 被 注释 到 公共 数据 库 , 这 些 
序列 可 能 是 一 些 非 编码 序列 或 者 新 发 现 的 基因 序 
列 ,也 可 能 是 因为 某 些 unigene 长 度 过 短 并 未 能 与 
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图 7 











应 用 邻接 法 基于 核 背 酸 序 列 构建 的 花椒 罕 吉 丁 雌雄 成 虫 触 角 转 录 组 中 5 个 AzanCSP 基因 





与 其 他 鞘翅 目 昆 虫 CSP 基因 的 系统 发 育 树 (1 000 次 重复 ) 


Fig. 7 Phylogenetic tree of five AzanCSP genes in the antennal transcriptome of female and male adults of 


Agrilus zanthoxylumi and CSP genes from other coleopteran insect species based 


on nucleotide sequences by using neighbor-joining method (1 000 replicates) 
CSP 基因 来 源 物种 Origin species of CSP genes; AzanCSP ; EAE T Agrilus zanthoxylumi ; AmalCSP ; 苹果 小 吉本 Agrilus mali; AplaCSP ; iL: 





YET Agrilus planipennis; TmolCSP ; 黄粉 虫 Tenebrio molitor; AchiOBP ; 星 天 牛 Anoplophora chinensis; TcasCSP: JR $ 





4X Xx Tribolium castaneum; 








RdomCSP : 48$ Rhyzopertha dominica; DarmCSP; 华山 松 大 小 过 Dendroctonus armandi; DponCSP; "P Kx il TA 7]N Dendroctonus ponderosae; 
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PmacCSP ; $ ti E 3 IH P Pyrrhalta maculicollis; PaenCSP; 榆 
CbowCSP ; 大 猿 叶 
金龟 Holotrichia parallela; AcorCSP ; 铜绿 丽 金 
helophoroides. 












































公共 数据 库 中 的 已 有 基因 序列 比 对 到 ,此 外 ,也 与 公 
共 数 据 库 中 现 有 的 辅 她 目 昆虫 吉 丁 科 的 参考 序列 的 
数量 不 够 丰富 有 关 。 本 研究 发 现 少数 仅 在 花椒 罕 吉 
触角 中 存在 的 气味 结合 蛋白 基因 ,这 些 特 异 表 达 
的 基因 可 能 与 花椒 罕 吉 丁 物种 生境 中 所 特有 的 化 学 
物质 有 关 ,与 它 的 强 隐蔽 性 和 对 寄主 专 一 性 也 紧密 
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Pyrrhalta aenescens; AglaCSP ; 光 肩 星 天 牛 Anoplophora glabripennis; 


Colaphellus bowringi ; MaltCSP ; WAR Monochamus alternatus; GdauCSP ; Yh Eut Galeruca daurica; HparCSP ; Wi EAR 
f& Anomala corpulenta; AgraCSP ; 棉铃 e 

















Anthonomus grandis; DhelCSP ; 1E S 23 P Dastarcus 





关联 。KEGG Jd JE X) IE BUE ST ff n 
unigenes 可 能 参与 或 涉及 的 代谢 途径 的 注释 结果 分 
析 表 明 ,花椒 罕 吉 丁 信和 号 感受 与 传递 的 主要 场所 是 
触角 。 前 期 我 们 已 经 对 雌雄 成 虫 触角 中 这 些 靶 标 基 
因 的 表达 差异 进行 了 初步 探讨 ( 巩 雪 芳 等 , 2019 ) , 
本 研究 中 4zan0BP3 在 触角 中 高 丰 度 表达 (图 8: 
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花椒 罕 吉 丁 肉 雄 成 虫 触角 转录 组 中 嗅觉 相关 基因 的 RPKM fH 


Fig. 8 RPKM values of olfaction-related genes in the antenna of female and male adults of Agrilus zanthoxylumi 


A: AzanOBP ; B; AzanCSP ; C; AzanOR ; D; AzanIR ; E; AzanSNMP.  AnF ; WEE fi ffg Antenna of female adult; AnM ; HÈ E fih 


adult. 


A) ,与 前 期 的 研究 结果 ( 巩 雪 芳 等 , 2020 ) 不 一 致 ， 
主要 是 因为 本 研究 对 象 是 花椒 窄 吉 丁 触角 (不 包含 
头 部 ) ,而 前 期 目的 基因 (4zan0BP3 ) 在 触角 中 的 表 
达 研 究 中 的 研究 对 象 是 花椒 罕 吉 丁 头 部 (包含 触 
f8) , 且 根据 系统 发 育 分 析 发 现 Azan0BP3 与 畏 翅 目 
昆虫 的 性 信息 素 结合 蛋白 聚 类 到 同一 分 支 ,初步 猜 
测 并 推断 其 很 可 能 是 与 花椒 窗 吉 丁 性 信息 素 相关 的 
气味 结合 蛋白 ( 巩 雪 芳 等 , 2020) ,所 以 在 瞧 雄 成 虫 
触角 中 的 表达 量 差异 是 显著 的 ,在 肉 雄 成 虫 头 部 的 
表达 量 差 异 是 不 显著 的 。 

筛选 到 7 个 候选 OBP 基因 , 相 比 于 其 他 鞘翅 目 
昆虫 的 触角 转录 组 研究 ,本 研究 所 获得 的 嗅 党 基因 
数量 较 偏 少 ,除了 受 NCBI 基因 数据 库 (https: / 
www. nobi. nih. gov/) 中 参考 序列 的 数量 的 影响 外 ， 
还 可 能 是 因为 化 学 感受 基因 中 的 气味 结合 蛋白 和 化 





















































Antenna of male 














学 感受 蛋白 两 个 家 族 进化 速率 较 高 ,这 些 基 因 进 化 
很 可 能 与 识别 特定 植物 化 合 物 有 关 , 本 人 研究 中 的 人 研 


究 对 象 花椒 窄 吉本 为 单 食性 昆虫 , 目前 发 现 寄主 植 
物 只 有 花椒 一 种 ,这 可 能 造成 其 OBP 和 CSP 基因 不 
如 其 他 昆虫 的 数量 多 。 发 现 花椒 罕 吉 丁 OBP 并 非 
都 是 典型 的 OBP( 含 有 6 个 保守 的 半 胱 氨 酸 )( 巩 雪 
2$ 55, 2019) ,这 7 个 OBP 基因 中 ,只 有 5 个 OBP 基 
因 具 有 全 长 开放 阅读 框 ; 4Azan0BP1，Azan0BP2， 
AzanOBP3 , AzanOBP4 和 AzanOBPS5 的 编码 蛋白 具 
^H 6 个 保守 的 半 胱 氨 酸 位 点 ,符合 气味 结合 蛋白 的 
村 E ( Cys-X15-39-Cys-X3-Cys-X21-44-Cys-X7-12- 
Cys-X8-Cys, X 代表 任意 一 种 氨基 酸 ) ,为 典型 的 
0BP;Azan0BP6 的 编码 蛋白 序列 中 总 共有 10 个 半 
胱 氨 酸 位 点 ,是 非典 型 的 气味 结合 蛋白 ( atypical 
OBPs) ;Azan0BP7 的 编码 蛋白 具有 4 个 保守 的 半 胱 
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氨 酸 位 点 ,初步 鉴定 为 DimerOBP 类 。 

利用 RPKM 值 对 基因 表达 量 的 评估 结果 表明 
(图 8) , AzanOBP1 和 Azan0BP2 在 上 肉 成 虫 触角 中 不 
表达 ( RPKM = 0) ,在 雄 成 虫 触角 中 微量 表达 ; 
AzanOBP3 在 雄 成 虫 触角 中 高 丰 度 表达 (RPKM > 
1000) , f WE y, E fh fA rp rp E de iA ( RPKM = 
205. 25) ,推测 其 在 触角 中 高 丰 度 表达 的 OBPs 功能 
与 雄 成 虫 识别 同类 异性 释放 的 信息 素 或 寄主 植物 释 
放 的 挥发 物 有 关 。 

本 研究 首次 构建 了 花椒 性 干 害 虫 花 椒 窗 吉本 的 
雌雄 成 虫 触 角 转 录 组 数据 库 ,并 成 功 筛选 得 到 了 部 
分 咒 觉 相关 基因 ,为 将 来 利用 嗅觉 骤 标 无 公害 防治 
农林 害虫 提供 了 重要 分 子 信息 ;基于 KEGG 数据 库 
把 从 已 经 完整 测序 的 花椒 窜 吉 丁 触角 转录 组 中 得 到 
的 基因 目录 与 更 高 级 别 的 细胞 、 物 种 和 生态 系统 水 
平 的 系统 功能 关联 起 来 。 研 究 结 果 为 鞘翅 目 其 他 近 
缘 物 种 昆虫 的 触角 转录 组 的 研究 提供 了 丰富 的 转录 
本 基础 数据 信息 ,对 畏 这 目 吉 丁 科 非 模式 昆虫 的 研 
究 有 重要 的 推动 意义 。 
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